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ABSTRACT 

In this paper we search for a signature of a large scale bulk flow by looking for fluc- 
tuations in the magnitudes of distant LRGs. We take a sample of LRGs from the 
Sloan Digital Sky Survey with redshifts of z > 0.08 over a contiguous area of sky 
Neighboring LRG magnitudes are averaged together to find the fluctuation in magni- 
tudes as a function of R.A.. The result is a fluctuation of a few percent in flux across 
roughly 100 degrees. The source of this fluctuation could be from a large dipole mo- 
tion with respect to the LRG sample, or a systematic in our treatment of the data 
set, or the data set itself. A dipole model is fitted to the observed fluctuation, and 
the three flow parameters, its direction and magnitude: ah, 5b, Vb are constrained. We 
find that the flow direction is consistent with the direction found by other authors, 
with ctb ~ 180, 6b ~ —50. The flow magnitude however was found to be anomalously 
large with Vb > 4000km/s. The LRG angular selection function cannot be sufficiently 
taken into account in our analysis with the available data, and may be the source of 
either the anomalous magnitude of the flow signal, or possibly the entire fluctuation. 
However, the fluctuation indicates a flow direction very close to those found using 
other data sets and analyses. Further investigation with upcoming data is required to 
confirm this detection. 

Key words: large-scale structure of universe - cosmological parameters - surveys - 
galaxies: kinematics and dynamics - galaxies: statistics 



1 INTRODUCTION 

Recent years have seen a rapid growth of cosmological ob- 
servations and analyses that has brought forth the era of 
. prec ision cosmology led by missions like the WMAP satel- 
lite (ISpergel et alj I2OO7I : iDunklev etafl 120091 : lLarson et ail 
2010), the SDSS collaboration (see at |http://www.sdss.orgp , 
DEEP (http://deep.berkeley.edu) and many others. As im- 
portant and interesting as these measurements are, prob- 
ing the large scale structure of the Universe using many 
diverse and independent measurements remains of utmost 
importance. As much as we have learned in the last decade, 
cosmolo gy and even param eter estimation is not a closed 
subject ( Bridle et al.l 12003 ). Not only do we not know the 
nature of dark matter, but even the density and ampli- 
tude of its initial perturbations (usually parametrized by 
fi m and og respectively) remain an open question. The de- 
termination of the matter content is especially pro blematic 
l| Fukugita fc Peebles! [20oi : lJuszkiewicz et al.li2010l ) . 
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Further, there are some indicatio ns that suggest that 
the ACDM cosmology with WMAP |Larson et all |2010| ) 
central parameters may be problematic. Some examples 
are: Large-scale anomalies fo und in the maps of tempera- 
ture anisotropies in the C MB (Sa rkar et al.ll2010l ; ICopi et al.l 
l20ld : lBenriett et al 1l201ll) ; re cent es t imates of the large scal e 
bulk flow bvlWatkins et alj (120091): iFeldman et al.l |2010t ): 
iMacaulav et all |201ll ): lMa et al.l (|201ll ) are inconsistent at 
the nea rly 3a level with ACD M predictions; a recent es- 
timate l|Lee fc Komatsul 120101 ) of the occurrence of high- 
velocity merging systems such as the Bullet Cluster is un- 
likely at a ~ 6<r level; large excess of power in the sta- 
tistical clustering of luminous re d galaxies (LRG) i n the 
photometric SDSS galaxy sample (|Thomas et al.ll201ll ): evi- 
dence of larger than expected cross correlati on between sam- 
ples of galaxies a nd lensing of the CMB l|Ho et al.l 120081 : 
iHirata et al.ll2008l ); brighter t han expected type la Super- 
novae (SNIa) at High Redshift (|Kowalski et al.ll2008l ); voids, 
especially smaller ones (~ 10 Mpc) are observed to be 
much emptier than predicted |Gottlober et al.ll2003l ): the 
predicted shallow low concentration and density profiles 
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of Cluster Haloes disagree with observatio ns which indi- 
cate denser high concentration cluster h aloes (|de Blokll2005l ; 
iGentile et alj2005h : lKovetz et all l|201Ch find a unique direc- 
tion in the CMB sky determined by anomalous mean tem- 
perature ring profiles, also centered about the direction of 
the flow detected above. 

The amplitude and growth of cosmological fluctua- 
tions on large scales are closely related to the CMB dipole, 
which reflects the bulk flow (BF) of the local group and 
can be used to test cosmological models. Large scale ve- 
locity i _surv^vs_liave beejiund^ group s 
(e.g. iGiovanelli et all lll997l); iHudson et alJ d 19991 l2004) ; 



iMasters et all <|2006MSpringob et al l ll2007h:lspringob et all 
d2009l)1 and r ecent analyses by IWatkins et al. I i|2009l ); 



iFeldman etafl (|2010h of the newest proper distance mea- 
surements show that virtually all velocity survey analyses 
show a consistent large-scale BF. It appears to have a mag- 
nitude J> 400 km/s on scales of 100/i _1 Mpc, which disagrees 
with ACDM WMAP predictions at the ~ 3er level. Analysis 
of the flow l|Feldman et al.ll2010T ) suggests that if it is due 
to a gravitational potential flow, then the sources (i.e. over- 
and under-densities) must be on scales > 300/i _1 Mpc. The 
direction of this flow is close to the Galactic disk, Galactic 
longitude ~ 295° and longitude ~ 10° with error of ~ 5°. 
Assuming that this is a potential flow, we expect to see 
over-densities in the flow direction and under-densities in 
the opposite direction. The responsible stru ctures must be 
far aw ay from us since the flow is very cold (|Feldman et al.l 
l2010l ). 

The motion detected in IFeldman et"all (|201Ch is not 
due to nearby sources, such as the Great Attractor (dis- 
tance of ~ 40/i _1 Mpc), but rather to sources at greater 
depths that have yet to be fully identified. The largest known 
mass concentration, the Shapley supercluster, does not seem 
to be massive enoug h to cause a flow of th i s mag nitude 
l|Ravchaudhurvl 1 19891 ) . Following iTullv et al l (|2008l ). it is 
more likely that the flow arises both from various mass con- 
centrations in the Galactic y-direction as well as under-dense 
regions in the opposite direction. Currently, there is no sur- 
vey in existence that is deep enough to resolve the source(s) 
of the flow. Nor is there any data that probes these scales 
(> 300ft _1 Mpc) to see whether there are any large mass 
concentrations and voids in these directions. 

iKashlinskv et al l l |200Sl . l201Ch found a larger amplitude 
BF, coined the Dark Flow, on scales out to 300 /t -1 Mpc 
by using the kinematic Sunyaev-Zeldovich effect (kSZ) in 
the CMB to measure the bulk velo ci ty of about 120 X-ra y 
clusters. Like IWatkins et"ai1 d2009h : IFeldman et"al] (|2010t ). 
these papers claim that a BF at these depths, as determined 
by their studies, is difficult or impossible to explain within 
the framework of standard ACDM model of cosmology. Al- 
though these results utilize completely different methods to 
measure the BF and find similar results, kSZ and velocity 
surveys are both subject to problems of very large noise and 
systematics in the data. In particular it should be noted 
that converting the dipole of the kSZ signal to a velocity is 
model dependent: assuming a different cluster radial profile 
can give a c ompletely different result. However, oth er re- 
cent results l|Davis et al.ll201ll ; iNusser fc Davisll201l[ ) using 
the SFIH — h catalog of Tully-Fisher galaxies and also recon- 
structing the cosmological large-scale flows using the 2MASS 
(Two Micron All Sky Survey) redshift survey (2MRS) found 
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Figure 1. The magnitude fluctuation as a function of R.A.. The 
LRG r band model magnitudes in bins of R.A. and redshift are 
averaged over all values of the declinations, the mean magnitude 
at each redshift is then subtracted. The fluctuation is in units of 
"maggies" , see Eqs. 1 1 121 Each panel shows the different spectro- 
scopic redshift bins: top panel z c = 0.12, 2nd panel z c = 0.20, 3rd 
panel z c = 0.28, bottom panel z c = 0.36. The solid lines are the 
Galactic extinction corrected magnitudes, and the dotted lines 
are the magnitudes with no extinction correction applied. 
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flows roughly the same direction as mentioned above but 
with magnitude consistent with the WMAP ACDM model. 

In this paper we approach the problem in a different 
manner and instead look for fluctuations of observed mag- 
nitudes across the sky as could be induced by a bulk motion 
centered on us. We choose objects that we assume are dis- 
tant enough to be considered at rest compared to the CMB. 
Ideally these objects should be standard candles, therefore 
the best astrophysical candidate would be SNIa. Unfortu- 
nately there are not enough observations of SNIa, and their 
distribution across the sky is far from ideal for this investiga- 
tion. We choose instead to use luminous red galaxies (LRGs) 
which are suitable because they are assumed to form a single 
population of galaxies, which assembled at high-redshift and 
have been pa ssively evolv i ng sin ce. This assumption was in- 
vestigated bv lWake et alj (|2006h who found that the LRG lu- 
minosity function (LF) is consistent with purely passive evo- 
lution and a merger-free history. Therefore to use the LRGs 
as a proxy for standard candles we average out the effect 
of the intrinsic variability of their luminosities by averag- 
ing tog ether the magnitud es of neighboring LRGs. A recent 
study (jNusser et al.l I2OI1T ) presented an alternate method 
for detecting cosmological bulk flows from redshift surveys, 
using the observed dimming or brightening of galaxies due 
to their peculiar motion, a similar technique to that which 
we present here. 

The paper is organized as follows, in Section [2] we 
present our data set and how we determine the observed 
magnitude fluctuation, in Section [3] we present the theoret- 
ical model for magnitude fluctuations induced by peculiar 
velocities. Section [4] describes our method of analyzing the 
data, with the results presented in Sections 14. 21 and 14. 31 and 
discussion and conclusions of these results in Section [5] 

Throughout this paper we assume a concordance cos- 
mology of: -ffo=70km/s, fl m = 0.3, J1a = 0.7 and w — —1 
unless otherwise stated. 



2 THE OBSERVED FLUCTUATION 

We use the SDS S spectroscopic and pho tom etric LRG sam- 
ples, d escribed in lEisenstein et all <|200lh and lCollister" et al.l 
(2007). We select a rectangular area of the SDSS survey 
bounded by: 125° < R.A. < 245° and 8° < Dec. < 59°, 
and covering the redshift range: 0.08 < z < 0.40 for the 
spectroscopic sample, and 0.4 < z < 0.7 for the photomet- 
ric sample. The area was chosen to be contiguous and to 
be about a degree or so from the survey edges to ensure 
a uniform as possible coverage. We use the r band model 
magnitudes and we further filter the catalog by requiring 
that all the galaxies have a signal-to-noise ratio > 2 and the 
Galactic extinction in the r band x r < 0.1. Redshifts are 
corrected to the local group frame. 

To obtain enough signal we must average over enough 
LRGs in order to beat down the effect of their intrin- 
sic variability in lu m inosity. The BF dire ction found by 
IWatkins et"afl (<2009l ): lFeldman et al.1 l|2010l i is towards R.A 
= 180°, Dec. = —52°. Since the R.A. direction lies within 
the survey area, whereas the SDSS LRG Dec. distribution 
are more than 60° away, we choose to average the LRG mag- 
nitudes across declinations. For each redshift and R.A. bin 
we find the average difference in flux units, or "maggies" 
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Figure 2. The "model dependent" magnitude fluctuation as a 
function of R.A.. The difference between LRG r band model mag- 
nitudes and the "theoretical" catalog in bins of R.A. and redshift 
is averaged, over all values of the declinations. The fluctuation 
is in units of "maggies", see Eos. I3I41 Each panel shows a dif- 
ferent spectroscopic redshift bin: top panel z c = 0.12, 2nd panel 
z c = 0.20, 3rd panel z c = 0.28, bottom panel z c = 0.36. The 
solid lines are the Galactic extinction corrected magnitudes, and 
the dotted lines are the magnitudes with no extinction correction 
applied. 
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Figure 3. The magnitude fluctuation as a function of R.A. for the photometric sample. Each row shows a different photometric redshift 
bin: top row z c = 0.475, bottom row z c = 0.625. The two columns correspond to the two different ways of computing the fluctuation: 
left column is the fluctuation described by Eq. [T](as in Fig. [T] for the spectroscopic sample), the right column is the "model dependent" 
magnitude fluctuation described by Eq. [3] (as in Fig. [2] for the spectroscopic sample). 



(a linear measure of flux defined as lO -0 ' 4 ™), between the 
magnitudes in that bin and the mean magnitude at any an- 
gular position within the same redshift bin: 

df = io- - 4 ^-^) (1) 
C = <df> t] (2) 

where i refers to the R.A. bin, and j to the redshift bin, df is 
the difference in LRG flux/magnitude in maggie units and 
5m is the mean value of df in the R.A. -redshift bin ij . Fig.[T] 
shows this magnitude fluctuation in each redshift bin for the 
spectroscopic sample; the dotted lines show the fluctuation 
when the magnitudes have not been corrected for Galactic 
extinction. The left hand side of Fig. [3] shows the same for 
the photometric sample. 

Computing the observed fluctuation in this way is not 
optimal because we do not have a full sky area with which 
to compute (m r )., and the result will also be sensitive to 
the redshift distribution of the LRG within each R.A. bin. 
We would also like to compute the observed magnitude fluc- 
tuation in such a way that it could be analysed in terms of 
the theoretical perturbation of observed magnitudes due to 
peculiar velocities. Therefore because of these two issues we 
also calculate the observed fluctuation in R.A., redshift bins 
in the following way: 



-0.4(r, 



df = TO 

<C = < df >ij 



-m r ) 



(3) 
(4) 



where again the observed fluctuation 5^ is in flux rather 
than magnitude units. Eq[T]and Eq. [3] differ in their use of 



fh r or (m r )j- fh r are the expected apparent magnitudes for 
each LRG in our catalog, assuming a cosmology, no flows 
(homogeneous universe) and that all LRGs belong to a pop- 
ulation with the same mean absolute magnitude and spec- 
tral energy distribution (SED), whereas {m r ) ■ are simply 
the average LRG magnitude (over all angular positions) in 
redshift bin j. 

Fig.[2]shows this fluctuation in the spectroscopic sample 
for each redshift bin, and with and without Galactic dust 
corrections applied to their magnitudes, and the right hand 
side of Fig. [3] shows the same for the photometric sample. 

The method we use to calculate m r , the expected ap- 
parent magnitudes in the absence of any flows for each LRG 
in our catalog, is outlined below: 



rn r {z) = 51og 10 d L (z) + 25 + (M r ) lrg + K r 



(•>) 



where rh r and dc denote the magnitude and luminosity dis- 
tance in a homogeneous universe respectively, {M r ) lrg is the 
average LRG absolute magnitude in the r band, K r is the 
fc-correction, z is the redshift in our filtered LRG catalog. 
dL(z) is ca lculated using the usual equations as given in 
iHogd (|l999l ). The fc-correction is found by assuming that all 
the LRGs have an SED with the same shape. It is calculated 
for each galaxy from: 



-2.51o gl 



1 Jd\\f x (\/(l + z))r(\) 
(1 + z) J d\\f x (\)r(X) 



(6) 



where f\(X) is the assu med SED of the LRGs , in this case 
it is the elliptical SED in I Coleman et all l|l980l ). and r(A) is 
the SDSS r filter. 
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For calculating the mean absolute magnitude {M r ) lrg 
we use the SDSS LR G luminosity function (LF) measured 
bv lWake et alJ l)2006l l corrected to redshift zero. To calculate 
{M r ) lrg we simply integrate the LF, 4>(M) multiplied by M r . 



(Mr) lrg = 



/ 4>(M)MdM 
J(j}(M)dM 



(7) 



From the above LF we find that we need to average over 
around 1000 galaxies in each bin so we can be confident 
that the average absolute magnitude of LRGs in each bin i, 

{M r ). EE (M r ) lrg . 

Finally this "theoretical homogeneous universe" LRG 
catalog is subtracted from the real LRG catalog, to get df 
as defined in Eq. [3] 

We now discuss the assumptions made in Eq. [S] when 
calculating the theoretical catalog. First of all we assume 
a fiducial set of cosmological parameters to calculate the 
luminosity distance di. The choice of parameters has a 
predictable effect on the resulting observed fluctuation 
5m, shifting the whole flucutation upwards or downwards 
slightly, i.e. changing the normalisation, but it cannot pro- 
duce an effect which varies with R.A.. Secondly we assume 
the value we calculate for (M r ) lrg is in fact the true average 
LRG absolute magnitude in the bin. As long as there are over 
500 galaxies in each bin (and we have arranged the binning 
such that this true) and taki ng into account th e errors on 
the luminosity function from Wa ke et al.l (|2006f ), the error 
on this quantity is less than 1%, smaller than the SDSS pho- 
tometric errors. There should be no reason to expect {M r ) lrg 
to vary with position, however it could evolve with redshift. 
Finally to com pute the k correction we must assume the el- 
liptical SED in IColeman et"afl (|l980h is truly representative 
of the actual LRG SED. Again we find no reason that this 
assumption, if wrong, would cause an effect where < rh r > 
varied with R.A.. In fact u sing a complete l y diffe rent SED, 
the spiral galaxy Sbc from IColeman et al.l (|l980l ). causes a 
negligible difference in the results because the difference in 
fc-correction values between these two SEDs are small. 

It seems that all the assumptions made in calculating 
this theoretical LRG catalog could cause the same system- 
atic effect: a change in the normalisation of the plots. This 
normalisation could be a function of redshift: the difference 
between dc calculated with different cosmological parame- 
ters increases as a function of z, and {M r ) lrg could evolve 
with z. Crucially, however, none of these systematics are 
expected to cause a change in (m r ) with angular position. 

In this paper we want to constrain the parameters of 
the BF: Vb, «{,, 5b- The observed fluctuation as defined by 
Eq. [3] is the one which is compatible with how the predicted 
magnitude fluctuations are modelled, as described in Sec.0 
Therefore to avoid issues with the "nuisance" systematics 
outlined above we divide out the normalisation of the fluc- 
tuations. We do this by dividing out the value: MAX(5„)/2 
+ MIN(<5^)/2 from the fluctuations shown in Fig. [2] 

Looking by eye at Fig. [T] and [2] it seems the first red- 
shift bin might be contaminated by LRGs which are not at 
rest with respect to the CMB, or indeed by objects from the 
Main SDSS sample. It is not unexpected that "normal" flows 
in the universe could still be affecting the observed magni- 
tudes at these redshifts. Therefore in Sec. 14. 21 we repeat our 
analysis using all redshift bins, removing the lowest redshift 



bin, and removing the two lowest redshift bins to reduce the 
possibility of either of these two effects. 

The two main candidates to explain at least part of the 
observed fluctuation presented in Figures [T] and [31 outside 
of the context of a dipole flow, are: reddening by our Galaxy, 
and the variation in the selection of SDSS LRGs. These 
can also be considered to be systematics. It is unlikely that 
Galactic reddening is a major contributor because i) there 
is so little difference in the observed fluctuation when us- 
ing extinction-corrected and non-extinction corrected mag- 
nitudes, ii) a similar fluctuation is seen in the g, i and z 
bands i.e. there is no fluctuation in the observed colors. To 
see the variation with extinction in the r band, compared to 
the average extinction, analogously to Eq. [T]we calculate: 



dx 
5" 



10 



-M('r-W,) 



< dx >i 



(8) 
(9) 



where x r is the extinction in the r band at the position of 
the LRG. This fluctuation is plotted at the top of Fig. U 
Each different line style corresponds to a different redshift 
bin. Because the Galactic extinction varies only with angular 
position, and of course not with the redshift of the LRGs, 
the pattern of mean extinction verses R.A. is the same for 
each redshift bin: small differences occur only due to the 
different distribution of LRG angular positions within each 
redshift bin. 

The LRG angular selection function is more difficult 
to quantify. To get a handle on what the angular selection 
function of SDSS LRGs might be we looked at th e NYU 
Value-Added Galaxy Catalog (|Blanton et al.l 120051 . NYU- 
VAGC). The NYU-VAGC is a collection of galaxy catalogs 
cross matched to SDSS designed for the study of galaxy 
formation, evolution, and large-scale structure. As such the 
catalog contains detailed information about the angular se- 
lection function. We choose the parameter f gotten from the 
large-scale structure subsample as the value we use to esti- 
mate the average LRG completeness as a function of posi- 
tion. The f gotten parameter describes the fraction of tar- 
gets in the SDSS Main sample which have successfully mea- 
sured redshifts. The bottom of Fig. [4] plots the mean value 
of fgotten in each of the R.A. bins divided by the mean 
value of fgotten over the whole survey area. Although the 
variation in fgotten is similar to the fluctuation we observe, 
its amplitude is significantly smaller. To incorporate the in- 
formation contained within the fgotten parameter into our 
analysis we would also need to know how it varied as a func- 
tion of flux, and not just angular position. 



3 THEORETICAL MODEL OF MAGNITUDE 
FLUCTUATIONS 

The perturbation to the luminosity distance is defined as: 
5dh = dh(z) — dh[z) % where dh(z) is the actual observed lu- 
minosity distance (i.e. given the universe is inhomogeneous) 
and dL{z) is the luminosity distance calculated at the same 
observed redshift z but assuming perfect homogeneity. A bar 
above quantities indicates that they have been calculated as- 
suming, or are defined by assuming, homogeneity. 

The perturbation to the luminosity distance (to first 
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Figure 4. Top: Average Galactic extinction in the r band in each 
redshift and R.A. bin, the different lines correspond to the dif- 
ferent redshift bins. Bottom: Average angular selection (f gotten) 
as a function of R.A.. The variation in f gotten seems to follow 
the variation in the Galactic extinction, which is not unexpected 
e.g. less galaxies are likely to be detected in regions of higher 
extinction. 



order) due to peculiar velocities is (|Hui fc Greenell2006l) 

e -f C(l + Z) 2 



Ml 
d L 



Hd L 



(10) 



where v<. is the peculiar velocity of the emitter, v„ is the 
peculiar velocity of the observer, f is the unit vector in the 
direction from the observer to the emitter, z is the redshift of 
the emitter, H and d,L are the unperturbed Hubble parame- 
ter and luminosity distance respec t ively, calculated from the 
usual equations as found in lHogd (| 1999T ) . 

The relation between the luminosity distance and ap- 
parent magnitude m is: 



m = 51og 10 d L +M = + M 

Then differentiating the above equation gives: 

5 5d L 



8m = 



In 10 d L 



(11) 



(12) 



We assume that the LRG's form a rest frame (v e =0) so we 
can look at our motion relative to them, our dipole flow v , 
the magnitude fluctuation is now: 



8m = 



5 v ■ f / c(l + z 



In 10 



Hd L 



(13) 



The value of v„ • f will depend on the angular position of 
the LRG (R.A. and Dec, labelled as a and 8 respectively). 
We model the values by assuming they are due only to a BF 



of magnitude Vb and in direction ab,Sb, so we can calculate 
Eq. Q3] as a function of position (a,5) and redshift z. 

The top of Fig. shows the change in iQ-° A5m w fth 
redshift for a constant value of v ■ f = 300km/s. Three dif- 
ferent cosmologies are plotted: ACDM in red; CDM universe 
(f2 m = 1) in green; Open universe (fi m = 0.3) in blue. The 
value of the fluctuation \Q-° ASm [ s no t a strong function of 
the cosmological parameters because here (in Eq. I13|l they 
only affect the "geometric" distance factors. The value of 
v ■ f (and hence the value of the fluctuation) would be pri- 
marily driven by the cosmological parameters affecting the 
growth of structure, in the absence of more exotic expla- 
nations. This is why the choice of cosmological parameters 
in Eq. [5] does not have any significant effect on the "model 
dependent" observed magnitude fluctuation. The bottom of 
Fig.[5]shows the expected magnitude fluctuation in our LRG 
sample if there is a dipole flow of magnitude Ub=300km/s in 
the direction of at, = 180 deg. and (5b=-50 deg. in a ACDM 
universe. 



3.1 The expected signal 

We expect the magnitude fluctuations generated by a dipole 
flow to follow a cosine form (in the magnitudes) across the 
sky, and we can model this cosine with the following func- 
tional form: 



Smag = A cos(a + B) + C 



(14) 



where a is the R.A. direction. Putting this together with 
Eq. [13] we can understand the effect of the flow parameters 
(a;,, 8b, Vb) on the resulting cosine shape of Smag vs R.A. 
in the following way: 

• A is the amplitude and is affected by Sb and Vb. 

• B is the phase and is affected only by ctb- 

• C is the normalisation and is affected by 5b and Vb, but 
also the systematics. 

There is a degeneracy between St and Vb which both af- 
fect the amplitude of the cosine. It also should be noted 
that two flows with a relation between their bulk angles as: 
otbi = 0*2 + 7T and 5bi = —8b2 are indistinguishable. Hence 
we only consider flows with 5b < 0. This qualitative expla- 
nation of the signal also illustrates why we can divide out 
the constant C: it is the only part affected by the systemat- 
ics, and information on the flow amplitude and declination 
direction is still contained within A. 



4 FITTING THE FLOW 

4.1 Errors on the measured 5 m 

Naturally, even in the absence of measurement errors and 
peculiar velocities, the LRG magnitudes will have some dis- 
tribution within the redshift bin. Then in the limit where 
the redshift interval of the bin is infinitesimally small, and 
all LRGs have exactly the same SED and negligible redden- 
ing, then this distribution should follow the LRG luminosity 
function. If the individual magnitudes are modelled as hav- 
ing a Gaussian measurement error distribution with a width 
of a, and if instead the LF was a 8 function at one magni- 
tude, then the observed distribution of magnitudes in the bin 
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Table 1. Results from the spectroscopic catalog 
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Figure 5. Top: Expected magnitude fluctuation amplitude as 
a function of redshift due to a dipole flow with value v ■ r = 
300km/s. Calculated from Eq. ll3l thcn transformed into maggies. 
Three different cosmologies are plotted: ACDM in red; CDM uni- 
verse (flm = 1) in green; Open universe (f2 m = 0.3) in blue. 
Bottom: Expected magnitude fluctuation for our LRG sample in 
each redshift bin due to a flow with a magnitude of v b = 300km/s 
in a ACDM universe. 



would be a Gaussian of width a around this value. There- 
fore we can approximate the a of the error distribution by 
looking at the difference between the observed distribution 
of magnitudes in each bin and the expected distribution: 



2 

a s m = 



(x-xy 



N 



(15) 



where X = m r — fh r : the observed magnitudes minus the 
expected magnitudes, X is the mean of this distribution, 
and N is the number of LRGs in the bin. 



4.2 x 2 fitting 

We calculate the following \ 2 statistic: 

8rn ^rn, (Vmag , C^b, $b) 



X 2 (vb,a b ,S b ) = ^2 



(16) 



where the sum is over the redshift bins i and the R.A. bins 
j; 8^ is the observed magnitude fluctuation in bin i,j (see 
Section. [2]) and 8^ is the predicted value of the magnitude 
fluctuation in bin ag m is the error on the observed mag- 
nitude fluctuation in bin i,j found from Eg. 1151 

To calculate the predicted fluctuation, 8^, we calculate 
the 8m N using Eq. [13] (given the trial flow parameters) for 



No. z bin 


a b 




St 


Vb 


2 


177.0j 


2.7 
1.4 


-50±g 


>7000 


3 


179.0j 


1.3 
1.7 


-65ti°s 


6000±™° 


4 


185.5j 


1.9 
3.0 


-35^° 





each of the N LRGs in our catalog, and then average them 
in the same way as the data: 



8 P m = < 10" 



(17) 



The probability as a function of the flow parameters is 



calculated from the x as follows: 



P(v b ,a b ,8 b ) oc exp ( - — 



(18) 



The resulting one-dimensional likelihoods of the flow param- 
eters (after marginalising over the other 2 parameters) are 
shown in Fig. [6] The 1-a constraints on each flow parameter 
(after marginalising over the other parameters) are shown 
in Table [T] for the spectroscopic catalog. 

Fig. [7| shows the magnitude fluctuation data in two 
highest redshift bins with the best-fit model from Table [1] 
over-plotted. The best-fit flow model was taken from the fit 
including all four redshift bins. 



4.3 Direction fitting 

We also attempted a separate approach to find the direction 
of the flow by fitting a vector of the directions of the LRG's 
weighted by each galaxy's magnitude (or maggie). We find 
the weighted mean R.A. of the fluctuations by weighing the 
R.A. of each LRG with the magnitude m™ or with its "mag- 
gie" Kr ' 4 ( m "- <m "» where is the r-band magnitude 
of each LRG indexed by n and < m" > is the mean of 
the distribution. We have done that for the whole distribu- 
tion (all redshifts) and separately by looking at the same 
redshift regions as described in Fig. [1] we also calculated 
the R.A. direction in various Dec. regions. The results are 
not very sensitive to the binning schemes and agree by and 
large with the analysis above. We found that the most likely 
R.A. is 178 ± 11° in good agreement with the binned analy- 
sis m^sented_liere and w it h the velocity analyse s presented 
in 



presented ncrc and wit n tne velocity analyse s presentea 
Watkins et all d2009h ; iFeldman et all (|2010l ): (Ma et all 



20111) ; iKashlinskv et al.1 ( |200^L l2010h ." 



5 DISCUSSION 

Our goal in this paper was to search for a signature of a bulk 
flow by looking for fluctuations in the magnitudes of distant 
LRGs. 

Figures [T] [5] and |3 show that we do find a coherent fluc- 
tuation in the LRG magnitudes which seems to be fairly 
independent of their redshift. The cosine-like shape of the 
fluctuation matches that which would be expected if this 
observed fluctuation is due to a dipole motion with respect 
to the LRG sample. The maximum of the fluctuation is also 
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R.A. (degrees) Dec. (degrees) velocity, (km/s) 

Figure 6. Onc-dimcnsional likelihoods of the flow parameters given by P a exp(— x 2 /2), after marginalising over the other 2 parameters 
for the spectroscopic catalog. The red dashed lines correspond to when all 4 redshift bins are used in the analysis; solid green lines when 
only 3 bins are used; solid blue when only 2 bins are used. 



located roughly at the R.A. direction of a bulk flow found by 
other authors (170±10°), however its amplitude is more than 
an order of magnitude larger than what would be expected 
from Eq. [13] after assuming a "reasonable" flow (Fig. [5]). We 
find a magnitude fluctuation on the order of a few percent 
in flux for the spectroscopic sample, and about a percent for 
the photometric sample. Comparatively one would expect a 
fluctuation of less than 1 percent if the flow had an ampli- 
tude of less than lOOOkm/s. Figures [TJ [2] and [3] do also show 
the amplitude of the observed fluctuations decreases with 
redshift, as expected from Eg. 1131 and illustrated by Fig. [5] 

Our results are extremely sensitive to systematic effects, 
which is an obvious alternate explanation for the magnitude 
fluctuations we find. Possible effects could include: photo- 
metric zeropoint variation, correlated errors in the Galactic 
extinction correction applied to the magnitudes, and selec- 
tion effects on the LRG sample. 

The SDSS imaging data is acquired in a continuous 
scan, and each scan is obtained along a stripe. The survey 
strategy means the data we average together end up coming 
from different stripes. However, it would be very unlikely 
that a zeropoint variation could produce such a coherent 
fluctuation after averaging across different stripes. 

Galactic extinction on the other hand will have an ef- 
fect which is correlated with the angular position of the line 
of sight. Figs. [T] and [2] show that there is very little effect on 
shape of the observed fluctuation when no Galactic extinc- 
tion correction is applied to the magnitudes. This implies 
that an error in the Galactic extinction estimate is very 
unlikely to be the sole source of the fluctuation. In Fig. [4] 
we plotted the average Galactic extinction in each bin. The 
shape of the extinction variation is extremely similar to our 
observed magnitude fluctuation, but its amplitude cannot 
fully explain the entire signal that we find. 

The LRG angular selection function, however, is likely 
to have a big effect on our results. The angular selection 
varies in a similar manner to the fluctuation we observe. It 
would seem from comparing the bottom panel of Fig. [4] to 
figures \T\ and O that this could only contribute to about 1 
percent of the observed fluctuation. It is not clear however 
how to take account of this in our analysis, one would need to 
know the angular selection as a function of LRG magnitude 
before precise constraints on the magnitude of the flow could 
be made. 

We fitted a flow model to the observed fluctuation, 



and constrained the three flow parameters, its direction and 
magnitude: oib, St, Vb- We found that the flow R.A. direction 
was consistent with the direction found by other authors. 
The R.A. direction of the flow was the best constrained pa- 
rameter, unsurprising because it is the only parameter where 
the rough best-fit value can be inferred from Fig.[T]and [2]by 
eye. It is also unsurprising, given the amplitude of our signal, 
that we find an anomalously large flow. We do not attempt 
to perform the fit using the photometric sample. Although 
it has the advantage over the spectroscopic sample in that 
it has over 6 times as many LRGs, and all of them are at 
redshifts where we are most confident that they are truly at 
rest with respect to the CMB, it is non-trivial to account 
for their redshift errors in the analysis. At this stage given 
the small size of the survey area and its direction we do not 
think this is a worthwhile exercise. 

We repeated the fit using the spectroscopic sample leav- 
ing out the lowest, and second lowest redshift bins, since 
these redshift bins may contain LRGs that are not at rest 
with respect to the CMB. As Fig. [S] shows these are also 
the bins which should contain the most signal. We find that 
when leaving out these bins some of the constraint on the 
declination direction of the flow is lost but the right ascen- 
sion results are robust. We also found that the magnitude 
weighted right ascension directional fit also agrees with our 
binned analysis well. 

Since our results indicate that our local group is mov- 
ing with respect to the SDSS LRG sample, we infer that 
the source of our local motion is from scales comparable to 
the LRG scale , z = 0(0.1 — 0.2). This result agrees with 
iFeldman et"al] (|2010T ) that found that the sources responsi- 
ble for the bulk flow are at an effective distance of J> 300 
/i _1 Mpc, i.e. well within the horizon, and contradicting the 
su ggestion of a coher e nt flow on m uch larger scales claimed 
by IKashlinskv et all (120081 |2010|) and modeled as a tilted 
Universe (IKashlinskv etHn1l200Sl ; iMa et alJlioIlT ) . A note of 
caution must be made here, the uncertainties and system- 
atic effects on the magnitude of the flow suggested here, 
precludes anything more than a hint of a subhorizon flow. 
We have presented results using only the r band data, using 
either the g or i band data instead does not significantly 
alter our conclusions. 

Iltoh et all (|2010T ) used galaxy catalogs (not just the 
LRGs) constructed from the SDSS Data Release 6 and 
looked at the variations in the pixelised number counts, 
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Figure 7. Examples of the 4-redshift bin best-fit model (bottom row of Table [T] blue dashed lines) compared to to the data (green solid 
lines). The left-hand side shows the 0.24 < z < 0.32 bin and the right-hand side shows the 0.32 < z < 0.40 bin. 



a similar approach to that presented in this paper. They 
found the probability of f3 (~ iWfc/c) was consistent with a 
zero bulk flow, but also was not inconsistent with flows with 
/3 ~ 0.01. In fact, the galaxy sub-samples closest to those 
used here (termed Northern Galactic Hemisphere - shallow 
samples) find values of f3 >~ 0.02, and are only marginally 
consistent with zero. This indica t es tha t the analysis both 
presented here and in lltoh et all |2010l ) may be seeing the 
same effect in the SDSS data: either the signature of a large 
scale flow or a systematic in th e data. However , in a sim- 
ilar method to lltoh et all (|2010j ). iBlake fc Wall (|2002l ) an- 
alyzed the radio galaxy distribution from the NRAO VLA 
Sky Survey and detected a cosmological dipole anisotropy 
that is consistent with the CMB dipole in both amplitude 
and direction. 

We also note that a non- negligible detection of a fluctu- 
ation in galaxy magnitudes should also show up as a detec- 
tion of excess l arge s cale power in the clustering of galaxies. 
I Thomas et alj (|201ll ) used a similar data set to the one em- 
ployed here, the MegaZ DR7 photometric LRG sample, and 
found excess power on large scales of between 2-4<r signifi- 
cance by looking at the angular clustering of the LRGs. 

iNusser et al.l (|201ll ) also used galaxy luminosities to 
study the bulk flow. Although they found a bulk flow with 
a similar direction to that found here, and by other authors, 
they do not find the flow to have an anomalously large mag- 
nitude as compared to that expected from ACDM. This con- 
tention between different measures of the bulk flow using 
kSZ, galaxy distances, galaxy luminosities shows that fur- 
ther studies of our local velocity field are vital to test the 
current paradigm. 

Upcoming surveys should be able to confirm or repu- 
diate this detection. The Large Synoptic Survey Telescope 
jLSST Science Collaborations et all 120091 . LSST) is po- 
tentially the best candidate for two reasons: its ability to 
discover of order 50,000 photometric SNIa per year, and its 
wide area coverage (20,000 sq.deg.). Most of these SNIa will 
not be able to be followed up spectroscopically, however 
it is expected that photometric redshifts of SNIa observed 
by LSST could reach a precision of just 1%. Ideally a 
multi-wavelength approach, which would necessarily involve 
different targets and different instruments, would provide 
the strongest evidence in confirming the existence of this 
signal. Further, the whole sky data from WISE (Wide-field 



Infrared Survey Explorer) satellite <|Wright et al.ll2010| ) will 
become available in the next few years. It will be sensitive 
to 10% overdensities out to ~ 700/i~ 1 Mpc (z ~ 0.2) and 
although it does not measure redshifts, it may be used to 
search for over- and under-densities in the direction and 
anti-direction of the flow. 
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